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Abstract 

We assume that all quark and lepton mass matrices have upper trian- 
gular form. Using all available experimental data on quark and lepton 
masses and mixing angles we make a fit in which we determine mass 
matrices elements. There are too many free parameters and our solu- 
tions are not unique. We look for solutions with small non diagonal 
mixing matrix elements. In order to reduce the number of free param- 
eters we assume that the matrix element (-M)i3 vanishes in all mass 
matrices. Such universal assumption was drown from considering dif- 
ferent numerical solutions. The lepton sector, due to large mixing 
angles and very small errors for charged lepton masses, is more re- 
strictive then quark sector. We present the solution in this case. The 
absolute values of neutrino masses are not fixed. The another pos- 
sibility of reducing number of free parameters was considered by us 
before. With the additional assumption motivated by SU(5) sym- 
metry which connects mixing in right handed down quarks with left 
handed charged leptons we get a solution in which observed Cabibbo- 
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Kobayashi-Maskawa mixing for quarks comes mainly from non diago- 
nal terms in up quark mass matrix. 
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The understanding of masses and mixing for quarks and leptons, after 
discovery that neutrinos are massive [1J, is one of the most interesting and 
hot problems [2] in the standard model. In this paper we shall consider rather 
unusual idea of mass matrices in triangular form. Such idea was proposed in 
[3]. In that paper, taking into account very specific additional assumptions, 
mass matrices for up and down quarks, neutrinos and charged leptons were 
given together with mixing matrices for left and right handed components of 
fields. We considered model which neglects CP violation as well as the one 
that took it into account in quark sector. Neutrinos, contrary to now popular 
assumption of being Majorana particles, were treated as Dirac particles. 

In this paper we also assume that the neutrino is a Dirac particle. We 
believe that because there is no positive result of double f3 decay experiment 
it is still a viable hypothesis. It could be worthwhile (smaller number of 
parameters) to study consequences of this assumption. We have shown in 
[3], first considering real matrices, that inclusion of CP violation for quarks 
was not a problem. In this paper we will consider real matrices neglecting 
CP violation because we believe after experiences from [3] that it will be not 
difficult to include it. In addition (M) 13 element of lepton mixing matrix 
seems to be very small (if at all different from zero) and we do not have any 
information about CP violating phase. 

How we can justify the assumption of using mass matrices in triangular 
form. We believe that left handed SU{2) doublets and not right handed 
singlets are responsible for mixing. After spontaneous symmetry breaking 
our basic objects that we have in lagrangian are mass matrices for up and 
down quarks, neutrinos and charged leptons. When we make a decomposition 
of mass matrix into diagonal matrix and two unitary matrices, taking into 
account that we have strongly hierarchical mass spectrum for up and down 
quarks as well for charged leptons and not so strongly for neutrinos, we will 
find out that matrix elements above diagonal influence mostly mixing for left 
handed components and below diagonal right handed components. If we want 
to have mixing determined by left handed components (active in interactions) 
we assume that in our mass matrices matrix elements responsible for mixing 
of right handed (below diagonal) components vanish. In this way we end up 
with upper triangular mass matrices that (with hierarchical mass spectrum) 
give very asymmetrical mixing much stronger for left handed then for right 
handed components. In this paper we continue to use the assumption of 
triangular form for quark and lepton mass matrices but we do not take into 
account all other specific assumptions we made in [3j. The masses of quarks 
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and leptons and the corresponding errors are very different spanning many 
orders of magnitude. To cope with this and the fact that for neutrinos we 
do not know masses but only mass squared differences we will use, as in [3] 
X 2 function (the least square method). Knowing measured mixing matrices 
for left components and their errors and masses of quarks and leptons with 
corresponding errors we will try to calculate mass matrices for quarks and 
leptons assumed in upper triangular form. 

In general we have too many free parameters and we can not find solution 
that is unique. We give examples of solutions with relatively small non 
diagonal matrix elements in mass matrices. The question is can we reduce 
the number of free parameters by assuming that some of them vanish. It is 
not enough to calculate formally number of free parameters. In the case of 
neutrinos and charged leptons it is not possible to get solution (with small 
X 2 value) when neutrino mass matrix is diagonal and whole mixing is in 
charged lepton sector in spite of the fact that formally we have enough free 
parameters to get a solution. We can get analogous solution for up and down 
quarks. We discuss in detail the different cases for quarks and leptons when 
we have enough free parameters and could not find the solution with small 
X 2 value. It seems that lepton sector is more restrictive then quark sector. 
The mixing for leptons is much stronger then for quarks and charged lepton 
masses are known with high accuracy (specially electron and muon masses) 
in comparison with up and down quark masses. As a result of our numerical 
analysis we have found out that when in addition to the assumption of upper 
triangular form of mass matrices we assume that in neutrino and charged 
leptons mass matrix matrix element (M)i3 vanishes we can find a solution. 
This assumption is generalized and is made for all multiplets of up and down 
quarks neutrinos and charged leptons so they are all treated in the same way. 
Such assumption can be justified by the fact that measured mixing matrices 
for quarks and leptons have matrix element (£7) 13 very small (smaller then 
the others, in the case of neutrinos it is not well known but very close to zero 
or even zero) so we expect that mixing of first and third generation comes as 
a result of mixing of other generations. Absolute values of neutrino masses 
can not be fixed in this way. 

At the end we will give rather exotic example of solution in which num- 
ber of free parameters is reduced by 3 by a relation between different mixing 
matrices. We thought at the beginning that it could be an alternative to 
the assumption of vanishing mass matrix elements. We will show that it 
is possible to get solutions for triangle mass matrices when we assume con- 
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nection between left and right handed mixing matrices. It is suggested by 
£77(5) symmetry condition that right handed components of down quarks 
transform in the same way as left handed components of charged leptons 
(belonging to the same £77(5) multiplet). The errors for electron and muon 
masses are very small so the mixing of charged leptons is rather small and 
that has a consequence that the contribution of mixing to CKM matrix from 
down quarks is not big and comes mostly from the mixing of up quarks. 
With our additional assumption for all mass matrices (M)i3 = we can not 
find such solution. 

Let us introduce the notation for quarks. We will follow the notation used 
in [1] and repeat introductory formulas. The £77(3) * SU(2) * U(l) gauge 
invariant Yukawa interactions for quarks are given by 

- L Y = QiiY^ijURjH^ + Qi(Yd)ijdRjH + h.c, (1) 

where Qi denote the ££7(2) doublets of left-handed quarks and ur, cIr are the 
right-handed up and down-type quarks respectively. The Yukawa couplings 
Y u and Y d are 3x3 matrices (i,j are the generation indices) and H is the 
£77(2) doublet Higgs field. After the electroweak symmetry breaking, these 
Yukawa interactions lead to the following quark mass terms 

- L m = uuiM^ijURj + d Li {M d ) i jdR j + h.c, (2) 
(M u )ij = (Y u )ijV, 

(M d )ij = (Y d )ijv, 

where v is vacuum expectation value of the neutral component of the Higgs 
field H. 

Mass matrices M u and M d (considered by us as basic objects) can be 
diagonalized using two unitary matrices U and V 

M u = U u M°Vl (3) 
M d = U d M°V d \ (4) 

The diagonal matrix elements (of and M®) correspond to the exper- 
imentally observed mass eigenvalues. The matrices U and V describe mixing 
between states with definite flavor u, d and with definite mass u', d! for left 
handed and right handed components 

u L = U u u' L , 
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u R = V u u' R , (5) 
d L = U d d' L , 
da = Vdd R . 

The generation mixing in the charged weak current after expressing in 
terms of fields with definite mass is described by the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix [I] which consists of two unitary matrices 

Uckm = UlU d . (6) 

In the similar way, assuming that neutrinos are traditional Dirac particles, 
we have mixing for the charged weak current in the lepton sector described 
by Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix [S] 

Umns = U}U U . (7) 
The relations for neutrino mass M v and charged leptons M\ are: 



M v = U v M°Vl (8) 
M l = UtMfV?. (9) 

As before diagonal matrix elements of M® and Mf 1 correspond to the ex- 
perimentally observed mass eigenvalues. 

We have got experimental data for CKM and PMNS matrices. We will 
use for CKM matrix and the errors of matrix elements values given in [6]. 
The values and the errors for up and down quark masses are taken from the 
[7J. The parameters of PMNS matrix and corresponding errors are taken 
from recent updating of the fit pj taking into account recent results from 
KamLAND and MINOS experiments given in [9]. Mass squared differences 
of neutrino masses are taken from the same fit [8j and masses of charged 
leptons and their errors are taken from [7] . The masses of charged leptons are 
known with extremely high accuracy in comparison with quark and neutrino 
ones. For neutrinos only mass squared differences are known. That was 
the reason why in [3], in order to fit mass matrices for quarks and leptons 
in triangular form to experimental data, we used numerical minimalization 
of x 2 function taking into account experimental errors. The corresponding 
modified function was used for PMNS matrix and mass values for neutrinos 
and charged leptons. Minimalization of \ 2 gives most probable value of 
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parameters or could give numerical solution to equations (as it was used 
in [3J. When there are too many free parameters we can get one of many 
possible solutions. In [3J we have also assumed that we have upper or lower 
triangular matrices for up, down, neutrinos and charged leptons and that 
matrices diagonalizing right handed components are the same in the weak 
isotpic spin multiplets. 

In the present paper we relax these additional assumptions. The left 
handed components of all quarks and leptons take part in charged weak cur- 
rent interactions. We want to have a situation that mixing in active left 
handed components determines the whole mixing. We assume that mass 
matrices for up, down quarks, neutrinos and charged leptons are upper tri- 
angular matrices. When we have upper triangular mass matrices with hi- 
erarchical diagonal elements decomposition into two unitary matrices and 
diagonal matrix we have very simple situation. Mixing matrix elements 
for left handed components U are determined to the first order by the 
ratios of non diagonal and corresponding diagonal mass matrix elements 



((U) 12 ~ (M) 12 /(M) 22 , ((E7) 23 ^ (M) 23 /(M) 33 , (([/)ia ^ (M) 13 /(M) 33 ) 



and matrix elements for right handed components V are determined by 
left handed mixing elements multiplied by corresponding ratios of diagonal 



masses ((M) n /(M) 22 , (M) n /(M) 33 , (M) 22 /(M) 33 ). Ortogonality of matri- 



ces U and V must of course be taken into account. When we have hierarchical 
masses the mixing in V is determined by U and much smaller. All that will 
clearly be seen in our example of the first solutions. In [5J we used artificial 
errors for electron and muon masses. We want to remind that masses of 
electron and muon are known with extremely high accuracy in comparison 
with all other particles. It seems that error in the mass of electron is compa- 
rable with the expected naively mass of heaviest neutrino. One can expect 
that these small errors can give some limitations on mixing in charged mass 
matrices because when we diagonalize them you have to end up with the 
values within small errors. We will try to see what are the consequences of 
small errors in masses of charged leptons in comparison with quarks. 
Following [3J for quarks we will minimalize the function 



X 
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(AU C K Mij ) 2 
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There are 6 free parameters (3 diagonal and 3 non diagonal) for upper 
triangular matrix and 6 parameters for upper triangular matrix M u (in 
principle we could also use lower triangular matrix for M u like in [3]). As 
we mentioned before the values m d * p = (md,m s ,mb), rrf^f = (m u ,m c ,m t ) 
and corresponding errors Am^ p , Amjjf are taken from the Particle Data 
Group [7] and U^km an d corresponding errors from [6] . The values of U C X ^ M 
are expressed in terms of angles to have unitarity satisfied to high degree 
and then we use calculated in this way matrix elements. For neutrinos and 
charged leptons we use the same procedure with obvious modifications. In 
the x 2 function we fit Uckm or Upmns matrix elements and masses of quarks 
or leptons using 12 free parameters ( 2 mass matrices, 3 diagonal and 3 non 
diagonal elements). There are more parameters we want to determine then 
independent fitted quantities so we can not uniquely determine them. Any- 
how, we do not understand the mass scales and big mass ratios for different 
generations of quarks. So it could happen that just because of big difference 
in mass scales non diagonal matrix elements are small in comparison with 
diagonal ones (for quantum systems very small mixing between very different 
energy levels) giving rather small mixing in Uckm matrix. We will look for 
example of this type of solution. 

We will start with Uckm matrix and up and down quark mass matrices 
and we consider the fit with relatively small parameters (elements of mass 
matrices are given in MeV) 
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The x 2 is practically zero (we get the number of order 10~""J so we can 
treat these parameters as a numerical solution of Eqs. (3,4) and Eq. (6). We 
see that relations between matrix elements of U and V mentioned before and 
mass matrix elements are satisfied. The obtained solution is not so different 
in character from that obtained with different additional assumptions in [3] . 
One can easily give (because we have 3 free parameters) examples with bigger 
non diagonal matrix elements giving the same Uckm mixing matrix. 

When we look for small mixing parameters in case of neutrinos and 
charged lepton mass matrices we get (in eV) 




0.82188 
-0.44863 
0.35107 

/ 0.00221 



0.00529 
0.01050 


0.56353 
0.55008 
-0.61632 



0.00327 
0.03329 
0.03502 

0.083388 
0.70437 
0.70491 



M, 



V,] 
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/ 0.98611 
0.16603 
V 0.00451 





0.00899 


-0.16461 
0.97333 
0.15975 






0.04904 

0.02214 
-0.15828 
0.98715 



(19) 



(20) 



(21) 



(22) 
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and (in MeV) 



I 0.510999 -0.03867 -0.07211 
105.658 -0.71477 

V 1776.99 



(23) 



/ 1 

0.00037 
V 0.00004 



-0.00037 -0.00004 \ 

1 -0.00040 
0.00040 1 / 



(24) 



/ 0.510999 




\ 

105.658 
1776.99 / 



(25) 
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1.77 x 10~ 6 
v -1.17 x 10~ 
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1.77 x 10~ 6 1.17 x 10~ 8 \ 

1 0.00002 
-0.00002 1 J 



(26) 



These are just examples of solutions with small non diagonal terms in 
mass matrices M u and Mi. We have relatively small mixing in U u and V u 
coming from M u and in U\ and V\ connected with M%. The mixing in Uckm 
and U pmns comes mainly from U d and U u connected with non diagonal terms 
of matrices M d and M v . For the quarks we have 3 free parameters and for 
leptons 4 (we do not know values of masses of neutrinos only mass squared 
differences) so we can find many very different solutions. In the CKM matrix 
Uckm = UlUd mixing can come from both matrices (the solution where most 
of the mixing in Uckm comes from U u is possible. That it is not the case for 
leptons. We do not have solutions where in Upmns mixing matrix most of 
the mixing comes from charged lepton mixing matrix U\. The mixing in the 
lepton case is much stronger then for quarks and quark masses are known 
with much less accuracy then masses of charged leptons. 

In the case of leptons situation is a bit different because of very small 
errors in masses of e and /i leptons. The diagonal masses are nearly equal to 
the measured charged lepton masses and non diagonal masses are really tiny. 
The biggest error is for r lepton mass. We have the strongest dependence on 
only one parameter namely element (M;)23 of charged lepton mass matrix. 
One can numerically try to find out how big this element could be still having 
solutions to our equations. We will give an example of the solution for leptons 
where we have relatively strong mixing because of "big" non diagonal matrix 
elements in mi mass matrix for charged leptons. 
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Mi 



I 0.00971 0.00278 -0.00038 

0.01493 0.02908 j (27) 
V 0.03942 

/ 0.51403 -11.3797 -9.29737 

105.647 -190.258 | (28) 
V 1766.71 



In this case U u and U\ are given by 



u, 



I 0.86365 0.50407 0.00323 

-0.39334 0.66988 0.62972 | (29) 

V 0.31527 -0.54513 0.77681 

/ 0.99411 -0.10826 -0.00519 
0.10708 0.98843 -0.10744 | (30) 

V 0.01676 0.10625 0.99420 

This is the solution that contrary to solution given in Eqs. (19-26) cor- 
responding to small non diagonal mass matrix elements gives relatively big 
non diagonal matrix elements for Mi. The mixing for charged leptons is still 
rather small in comparison with that for neutrinos. That makes that situation 
for leptons is different from that for quarks. In the case where non diagonal 
terms in charged lepton mass are relatively big the mixing in charged lepton 
mass matrix Ui is small and is only a small correction to U v being mainly 
responsible for the value of Upmns- Comparing with Eq. (27) with Eq. (19) 
we see that there is some increase in diagonal matrix elements in M u and that 
corresponds to bigger neutrino masses ni\ v = 9.1 meV, m2 V = 12.6 meV, 
m$ v = 49.8 meV in comparison with ni\ v = 2.2 meV, mi v = 9.0 meV, 
m$ v = 49.0 meV in Eq. (21). Contrary to quarks, in the lepton case we have 
rather small range of non diagonal terms in charged lepton mass matrix. The 
mixing given by Ui and VJ is small and the values of Upmns are dominated 
by non diagonal terms in M v and mixing matrix U v . 

Up to now we have considered examples of solutions with too many free 
parameters which we want to determine. One can restrict the freedom of 
possible solutions by some additional conditions e.g. as in [3] when the 
equality of V matrices in the same weak isospin multiplet was assumed. We 
will present mass matrices when we limit the number of parameters in mass 
matrices by putting 3 of them equal to zero. We will start with extreme cases 
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when when all vanishing mass matrix elements are in up or down quark mass 
matrix. In the case of quarks we get for diagonal M u : 



/ 5.13392 21.5375 16.0237 \ 
M d = 92.6034 175.388 
V 4196.3 / 



(31) 



and U d — Uqkm, whereas for diagonal M d we have: 



( 2.31044 -284.247 996.199 \ 

1218.37 -7167.79 
V 172348 / 



(32) 



and of course C/t = U C km- 

These are extreme cases when mixing in Uqkm matrix either comes from 
non diagonal terms in down quark mass matrix or from non diagonal terms 
in up quark matrix. We can also numerically study solutions when three zero 
matrix elements are distributed among two such matrices. It is easy to check 
that in spite of the fact that formally we have enough free parameters (3 zero 
matrix elements in M u and M d matrices) such solutions do not exist (x 2 is 
not small). When (M d )i 3 = 0, (M u )i2 = 0, (M u )\ 3 = or we change indexes 
u and d namely (M u )i 3 = 0, (M d )i2 = 0, (M,j)i 3 = there is no solutions 
(X 2 > 700 or x 2 > 1400). There is also no solution when (M d ) 23 = 0, 
(M u ) 13 = 0, (M M ) 2 3 = 0. That means that when we have number of free 
parameters equal to the number of experimental data we not always can find 
a solution. On the other hand when (M d )i 3 = 0, (M u )i 3 = and (M u ) 23 = 
or (M,j)i3 = 0, (M u ) 12 = 0, (M u ) 2 3 = we can easily find a solution. 

In the case of leptons when charged lepton mass matrix is diagonal (non 
diagonal elements are put equal to zero) we get: 



and Upmns = U v . It is also possible to have solution for M v with m\ v close 
to zero. 
We get 




0.00385 0.00510 0.00328 \ 
0.01100 0.03323 
0.03509 J 



(33) 



M v = 



( 10~ 12 0.00549 0.00329 \ 

0.01005 0.03337 
V 0.03494 J 



(34) 
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We want to stress that with diagonal mass matrix for charged leptons we 
still have one free parameter (we do not know absolute values of neutrino 
masses) and the smallest neutrino mass equal to zero is not excluded. 

In the case of leptons contrary to quarks, where we could have solu- 
tion with diagonal Md (Eq. (32)) but it is not possible (at least with high 
probability \ 2 > 170) to have a satisfactory solution that neutrino mass ma- 
trix is diagonal and the whole strong mixing that we have in Upmns comes 
from charged leptons. It is understandable from the discussion of mixing 
in charged lepton sector. With relatively small mixing in the charged lep- 
ton sector we can not reproduce very strong mixing observed in Upmns- 
The question is can we split three zeros between neutrino mass matrix and 
charged lepton mass matrix to have U u and U\ mixing matrices that repro- 
duce strong mixing in Upmns- To reproduce Upmns with two big mixing 
angles we need nonzero matrix elements (M v )i2 and (M v )23- Having only 
these matrix non diagonal matrix elements different from zero with all the 
other non diagonal matrix elements in M/ equal to zero we get \ 2 bigger then 
2. With (M v ) n and (M u ) 23 different from zero and (M,) 12 or (M,) 13 or (M,) 2 3 
we always get \ 2 bigger then 2. With 3 additional zeros in mass matrices we 
still have 1 free parameter connected with lack of scale for neutrino masses. 
Calculating formally number of parameters we should have a solution. The 
next possibility is to add one more non zero non diagonal matrix element. 
If we want to treat neutrinos and charged leptons in the same way then the 
choice is obvious. We assume that in charged lepton mass matrix non diag- 
onal matrix element (Mi) 23 is different from zero. In this case we can find 
a solution. That means that guided by numerical calculations we assumed 
that matrix elements (M u )i3 and (-Mi) 13 are equal to zero. We will also as- 
sume that for mass matrices of up and down quarks. It is the assumption 
corresponding for triangular matrices to that considered and advertised for 
quarks by Fritzsch, idea of having mixing only between neighbored flavors 
|10j . On the other hand from experiment we know that (Uc km) 13 is small 
and in {Upmns)i3 is small (not well known) and still consistent with zero so 
assuming that in basic mass matrices (M) 13 vanishes is not so unnatural. In 
some sense the situation is very simple. We ask is it possible to reduce num- 
ber of solutions by assuming that some mass matrix elements are equal to 
zero. If we want to add in a universal way (in the same way for up and down 
quarks, neutrinos and charged leptons) condition of vanishing mass matrix 
elements we do not get a solution with 2 zeros in every mass matrix. The 
next possibility is to consider one vanishing mass matrix element in every 
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mass matrix. Discussion of numerical solutions in lepton sector has helped 
us to localize these vanishing mass matrix elements. We will present a solu- 
tion with two additional zeros when matrix elements (M u )i3 and (Mi) 13 are 
equal to zero in neutrino and charged lepton mass matrices. In this case we 
get: 



M,, 




0.00381 
0.01125 0.02987 
0.03845 

0.8632 0.5047 0.01141 
-0.39654 0.6639 0.63403 
0.31242 -0.55183 0.77323 



/ 0.00455 







0.00983 




and (in MeV) 



M l = 



Ui = 



( 0.96535 
0.26096 
V 0.00118 



/ 0.51368 





-0.25834 
0.95498 
0.14584 



-10.7738 
105.655 




M, 



D 



f 0.99479 
0.10144 
V 0.01035 

0.510999 





-0.10197 
0.98961 
0.10133 



105.658 




\ 


0.0492 J 

0.03693 \ 
-0.14109 
0.98931 j 





-180.361 
1767.78 

0.00004 
-0.10186 
0.9948 




1776.99 



(35) 



(36) 



(37) 



(38) 



(39) 



(40) 



(41) 



1 0.0005 2.99 x 10~ 6 \ 

-0.0005 0.99998 0.00606 j (42) 
\ 1.06 x 10- 8 -0.00606 0.99998 / 

we have got a solution in which mixing matrix elements with value 0.1 in 
charged lepton mixing matrix U\ give substantial contribution to lepton mix- 
ing matrix Upmns- They are smaller then neutrino mixing matrix elements 
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but not as tiny as in examples of solutions given in the beginning (Eq. (24). 
This solution is not unique. It is possible to get solution (with accuracy 
10~ 13 ) with higher values of neutrino masses with slightly modified neutrino 
mass matrix and no change in lepton mass matrix parameters. One should 
remember that in the case of neutrinos we have one additional free parame- 
ter compared to quarks because we do not know absolute values of neutrino 
masses only mass squared differences. It is not possible to fix absolute val- 
ues of neutrino masses in this fit. We can also find a solution with smallest 
neutrino mass equal to zero with x 2 — 0.24. The value of \ 2 is higher then 
for solutions mentioned before, it can not be excluded, but is highly unsatis- 
factory. 

With the additional universal assumption (M) 13 = we can also look for 
solutions of this type for quarks. As we know from the previous discussion 
there is no problem with getting this type of solution. We know that it is 
possible to give a solution for quarks in which in addition to (M u )i 3 = also 
(M u ) 2 3 = 0. In this case we get: 




5.04798 13.0481 
94.18 





(43) 



U d 



\ 



0.99047 
-0.13761 
0.00578 



0.13773 
0.9896 
-0.04155 



2.92 x 1(T 6 \ 
0.04195 
0.99912 / 
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/ 5 \ 

95 
V 4200 J 



(45) 
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0.99997 
0.00732 
3.51 x 10' 



-0.00732 6.88 x 10~ 6 \ 
0.99997 -0.00094 
0.00094 1 / 



(46) 



and 



M u = 



I 2.25937 -113.728 
1244.82 

V o o 
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172500 / 



(47) 
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/ 0.99585 -0.09098 
0.09098 0.99585 
V 






1 



(48) 



M: 



u 



D 



( 



2.25 





\ 
1250 
172500 / 



(49) 



/ 



1 

0.00016 




0.00016 
1 
1 







(50) 



V 



We see from this solution that mixing in {Uckm)i2 is split, both mixing 
matrices Ud and U u contribute to Uckm- It is similar to neutrino and charged 
lepton case in Eqs. (36, 40). Matrix element {Upmns)v$ is reproduced by 
mixing in (M u ) 12 and (M^) 2 3 in mass matrices for up and down quarks. 
Unlike in the case of leptons mass matrix for up quark is not uniquely deter- 
mined. It is possible to give examples of solutions in which matrix element 
iMu)2z is different from zero ((M d ) 2 3 also changes) but mixing in (U u )%2 (that 
is bigger) is not strongly influenced. We have much stronger mass hierarchy 
for up and down quarks (with the mass errors relatively much bigger then in 
the case of charged leptons) and small mixing in Uckm in comparison with 
Upmns so the lepton sector is somehow more restrictive in our numerical 
analysis. 

We also want to give an example of solution in which number of free 
parameters is reduced by 3 (but still we have formally 4 free parameters) by 
the relation between different mixing matrices. We thought at the beginning 
that it could be an alternative to the assumption of vanishing mass matrix 
elements. On the other hand from the grand unification SU(5) symmetry 
(5 representation) we have a hint about connection between right handed 
components of down quarks and left handed components of leptons. This 
components being in the same SU(5) multiplet when there is a mixing be- 
tween families and in some way SU(5) structure is not completely lost in 
low energies [TTJ we would have the same mixing matrix for right handed 
down quarks and left handed charged leptons. We will now look what are 
the consequences of assumption Vd = U (see also [12]). It is not clear how 
well this relation could be satisfied because of corrections so we will assume 
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artificial error 0.001 (we know that mixing in C/j is rather weak). Minimaliz- 
ing x 2 function being sum of all terms corresponding to quarks and leptons 
and taking deviations of Vd from Ui with artificial error we can find solution 



/ 4.95067 5.4794 9.94122 \ 
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We get for this solution x 2 = 0.0086 with formally 4 free parameters. 
This solution is not unexpected. Small mixing in \J\ enforces by the condition 
Vd = Ui small mixing in V d and it means that mixing in Uqkm comes mainly 
from the mass matrix for up quarks. The obtained solution corresponds to 
the solution for M u and M u matrices given in Eq. (32) and Eq. (33). The 
mass matrices given in Eq. (32) and Eq. (33) with diagonal matrices for down 
quarks and leptons satisfy SU (5) condition as unit matrices. It is difficult 
to find a solution with a small x 2 value corresponding to the situation where 
most of the mixing in Uckm comes from down quark and only small part 
from M u mass matrix. It is not possible to find such solution with additional 
assumption that in all mass matrices matrix element (M)i 3 vanishes. 
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We have tried to determine mass matrices for quarks and leptons assum- 
ing that they have triangular form. For all particles it is assumed that mass 
matrices are upper triangular matrices. It means that in the decomposition 
of given mass matrix into diagonal and two unitary matrices left handed 
components of fields determine the mixing matrices. For hierarchical mass 
spectrum mixing in left handed components is much stronger then in the 
right handed components. As an input we use the values of matrix elements 
for CKM and PMNS matrices, masses of quarks and charged leptons with 
known errors and neutrinos squared mass differences with corresponding er- 
rors. With these assumptions one can not determine the elements of mass 
matrix uniquely using fit to experimental data. We have discussed examples 
of solutions for quarks and leptons with relatively small non diagonal matrix 
elements. It was stressed that that the masses of charged leptons are known 
with extremely high accuracy (the error in the determination of the mass of 
electron is comparable with naively expected highest neutrino mass). To- 
gether with very strong mixing in lepton sector, lepton sector seems to be 
more restrictive then quark sector. Considering various numerical solutions 
for leptons we came to the conclusion that in order to restrict the number of 
possible solutions and treat neutrinos and charged leptons in the same way 
we have to assume dynamical condition (M^)i 3 = and (Mj)i 3 = 0. We 
have given the solutions under these assumptions. Unfortunately, that does 
not fix the absolute scale for neutrino masses. Then we have extended this 
assumption to up and down quarks. The experimental justification for this 
assumption is that in Uckm and Upmns mixing matrices matrix element 
(U)i3 is much smaller then the other mixing matrix elements. We have also 
at the beginning as alternative to vanishing mass matrix elements considered 
other possibilities. With the additional assumption which connects mixing 
matrices for right handed down quarks and left handed charged leptons in 
the way suggested by SU(5) symmetry it is possible to find a solution for 
mass matrices with weak mixing for down quarks and charged leptons. The 
mixing is stronger for neutrinos and up quarks (mixing of left handed up 
quarks is mainly responsible for mixing in CKM matrix). With our addi- 
tional assumption for all mass matrices (M) 13 = we can not find such 
solution. 
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